General boundary conditions for the envelope function in multiband k ■ p model. 



A. V. RodinaU 

/. Physics Institute, Justus Liebig University of Giessen, 35392 Giessen, Germany 

A. Yu. Alekseev 

Institute of Theoretical Physics, Uppsala University, S-75108, Uppsala, Sweden 

■ Department of Mathematics, University of Geneva, 1211 Geneva, Switzerland 

(N . 

Al. L. Efros and M. Rosen 
^ 1 Naval Research Laboratory, Washington DC 20375, USA 

^ . B. K. Meyer 

0^ ' Physics Institute, Justus Liebig University of Giessen, 35392 Giessen, Germany 

We have derived general boundary conditions (BC) for the multiband envelope functions (which 
' do not contain spurious solutions) in semiconductor heterostructures with abrupt heterointerfaces. 

These BC require the conservation of the probability flux density normal to the interface and guar- 
antee that the multiband Hamiltonian be self-adjoint. The BC are energy independent and are 
characteristic properties of the interface. Calculations have been performed of the effect of the 
general BC on the electron energy levels in a potential well with infinite potential barriers using a 
coupled two band model. The connection with other approaches to determining BC for the envelope 
function and to the spurious solution problem in the multiband k ■ p model are discussed. 
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Q ; I. INTRODUCTION 

It is impossible to overestimate the role of the multiband effective mass approximation (MEMA) , which consists of 
the multiband k ■ p method together with the envelope function approximation (EFA), in the simulation of electronic 
0^ ' and optical devices formed from various semiconductor heterostructures, such as those with Type-I and Type II 
heterointerfaces, quantum well structures and superlattices, and one and zero dimensional semiconductor structures 
such as quantum wires and quantum dots. The calculations of energy bands and wave functions in such structures 
have often been considered as only being a matter of "energy band engineering" or "wave function engineering". 
However, although such calculations successfully describe many heterostructures (see, for example, Ref. [j]), a general 
\ theoretical description of structures with abrupt interfaces has yet to be developed. 

■ The MEMA gives an accurate description of the energy band structure of bulk semiconductors near the extrema of 
the first Brillouin zone. In this approach the carrier wave functions are expanded in the band edge Bloch functions near 
an extremum with "envelope function" coefficients, or components, that are required to vary slowly over a distance on 
the scale of the unit cell ao. The bulk MEMA Hamiltonian contains such parameters as band edge energies, interband 
momentum matrix elements and carrier effective masses. These material parameters are determined by fitting to 
independent experimental data or are obtained from first principles calculations. The MEMA is certain ^ appropriate 

\ for calculating the spectrum and carrier wave functions in the presence of any slowly varying potentialcl. r, p. 

However, abrupt heterointerfaces in semiconductor heterostructures break the periodicity of the crystal potentiallm. 
The nonperiodic part of the potential varies rapidly in an interface region of size a w oq. This mayJead to such interface 
effects as the admixture of remote band statesa, coupling of states with different symmetry^ - B, and the formation 
of surface localized statesa. Moreover, pEUoch functions in dissimilar materials can be significantly differentEl and do 
not even exist in the boundary regionstj. In light of the great success of the MEMA, it is important to develop an 
adequate description of semiconductor hetecostructures with abrupt boundaries using the MEMA to describe each 
material in the heterostructure. Laikhtmariij has shown that this is possible for heterostructures with a characteristic 
length L>a, provided that the boundaries are treated appropriately. The problem is then focused on the appropriate 
choice of boundary conditions (BC) for piece-wise defined envelope functions in these heterostructures. 

Conservation of the current normal to the interface is an obvious requirement dictated by conservation of probability 
density. In the EFA this condition reduces to the continuity across the interface of the normal component of the 
probability flux density J r averaged over the unit cellaO. In a heterostructure formed of semiconductors with 
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parabolic bands, which can be described by a single band effective mass approximation (EMA) (the simplest MEMA), 
this condition can be expressed as the continuity across the interface of: 

j T = -i£-(rf-f*f), a) 

2m V / 

where f is a unit vector directed normal to the interface, m is the effective mass, / is the single band envelope function 
and the prime denotes the derivative along r. The continuity of Eq. ([l]) is satisfied if 

fr = fl, ^ = A, (2) 

m r mi 

where the indices r and I refer to a point just to the right and just to the left, respectively, of the interface. These BC, 
expressing the separate continuity of the envelope function and the normal envelope velocity, / /m, first suggested in 
Refs. ia*e referred to as "conventional" BC. However, the conventional BC of Eq. (M^isweef ^cognized almost 



20 years ag to be a special case of more general BC that conserve J T across the interfao 




J! )=T tr ( J l ) , det [Ttr] = intffl - t 12 t 2 i = 1 (3) 
fr/m r J \fl/ m lj 

where the elements, Uj, of the,tcansfer matrix, Tt r , are characteristic of the interface0i, and, as seen from a more 
general derivation of Eq. (||)EjO, based on the requirement for the EMA Hamiltonian to be self-adjointlia, do not 
depend on the particular states and arp-|energy independent. For example, first principle calculations performed 
for GaAs / Al x Ga\- x As heterostructurespa show that the off diagonal elements tvi (which was shown to be always 
proportional to the small parameter a/l£3) and t^i are negligible for this interface. In this case the BC can be written 
in the form 

m<* r f r = mff h -A^^Jj^, (4) 

similar to the general BC first suggested in Ref. |l9L In general the parameters a and the off diagonal element t 2 i 
(which may be non-zero for other heterostructureiEl) strongly effect the penetration of the wave function into the 
barrier and the energy spectra in quantum wells (QWs). For example, this has been demonstrated in Ref. |2l] by 
analyzing the experimental data for the exciton optical spectra in narrow GaAs / AlGaAs and InGaAs/InP QWs 
using the EMA with modified BC similar to those given by Eq. (Q). However, the parabolic band approximation 
has a very narrow range of applicability and appropriate BC for the MEMA need to be developed for describing 
experimental data. 

For the multiband analysis of heterostructures, a natural extension of the conventional BC of Eq. (|^) is the continuity 
of each component of-theJV-component envelope function, 'J', and of each component of the normal envelope velocity, 
V T , at each interface .E3~EJ Taken together, these BC preserve the normal component of the envelope flux depsity, J T , 
averaged over a unit cell. The most general derivation of these BC, suggested by Barraff and Gershoni,EJ is based 
on the requirement for the MEMA Hamiltonian to be selt^idjoint. A general form of the BC within the MEMA for 
the planar heterostructures was suggested by Kisin et alS3. Those authors start from the conservation of J T at the 
interface, and introduce a transfer matrix connecting the components of the multiband envelope function, 'J, and its 
normal derivative, VP , on each side of the heterointerface. The general consideration of the 8x8 band model based 
on the symmetrical properties of the planar heterointerfaceEil has allowed us to gain important information about 
the elements of the transfer matrices. However, the BC introduced in Ref. [ll] may lead to an energy dependent 
connection between the components of \E r and V T on the left and right sides of the heterointerface, and, therefore, do 
not generally satisfy the requirement that the MEMA Hamiltonian be self-adjoint. 

The MEMA that take band coupling into account have a significant theoretical problem when describing various 
semiconductor heterostructures. The full expansion of the A-band envelope function in the N independent plane 
wave solutions of the bulk Hamiltonian corresponding to the same energy e includes solutions with a very-large wave 
vector lying outside the first Brillouin zoneEj'E£l. Such "wing" (evanescent )Ej or "spurious" (oscillatory)t£l solutions 
must be rejected in the perturbation k • p theory; in the bulk they are rejected as beimpjpchysical solutionsLZrEl 
However, their elimination from the multiband envelope function for heterostructuresE3Lirc3 : El leads to additional 
complications in the analysis of the BC. For example, the continuity of all N components of this truncated envelope 
function can_aot be satisfied if the interband matrix elements of the momentum operator p are not continuous across 
the interfaced. The same problem remains for the alternative approach, based on modifying the MEMA Hamiltonian 
by discarding those terms responsible for the large k solutionsOOEBEa. The use of different BC for the conduction 
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and valence band states (see, for example, Ref. [32]) automatically results in the dependence of the BC on e, and is 
justified only if band coupling in the MEMA Hamiltonian and bands nonparabolicity can be neglected. 

In our paper we consider the boundary condition problem within the MEMA for heterostructures with abrupt 
boundaries. We show that the most general requirement for the MEMA Hamiltonian to be self-adjoint implies 
continuity at the interface of the normal component of the envelope flux density vector J"' 3 . Here the indices a 
and (3 denote, in general, two arbitrarily chosen envelope functions \I/ Q and ^p, respectively. This allows us to 
formulate general BC using a new energy and state independent transfer matrix connecting the components of the 
envelope functions and their normal velocities on each side of the interface. It also allows us to find a general way 
of eliminating the large-fc solutions because these solutions make only a small (proportional to a/L) contribution to 
the envelope flux density. This procedure does not depend on the particular nature of the large-fc solutions (which 
we shall hereafter call spurious) and can be applied to boundaries with either finite or infinite potential barriers. 
The elements of these transfer matrices can be considered as characteristic surface or interface parameters and are 
additional to the bulk energy band parameters in the multiband k • p approximation. 

The paper is organized as follows. In Sec. |l| we derive the general form of BC for the multiband envelope function, 
that follows from the condition that the MEMA Hamiltonian be self-adjoint. In Sec. HI we derive a truncated form 
of general BC for slowly varying envelope functions which do not allo w sp urious solutions. To clarify the general 
procedure, we first consider, in detail, the two band Kane model. Sec. |v| gives analytical and numerical examples 
of the effect of the general BC on the quantum size energy levels in dots with infinite potential barriers. In Sec. |v| 
we discuss our results and compare them with those of other approaches to the BC problem within the MEMA and 
other k • p theories for heterostructures . 



II. GENERAL FORM OF MULTIBAND BOUNDARY CONDITIONS IN HETEROSTRUCTURES WITH 

ABRUPT BOUNDARIES 

Let us consider a semiconductor heterostructure made of M arbitrary shaped semiconductor layers with a charac- 
teristic length L. Single-electron wave functions, ^(r), in these structures are solutions of the Schrodinger equation 
with the microscopic Hamiltonian H m ia containing a crystal potential. This potential is "periodic" inside each of the 
bulk-like semiconductor layers, but it varies rapidly in the boundary regions with width a ~ ao. Assuming that a/L 
is small, we neglect the exact behavior of <f>(r) in the boundary regions and expand $ J ' e $(r € Rj) in the bulk-like 
interior region Rj of the j-th layer in the Bloch functions, u J nk , at a critical point, koj, of the bulk energy band 
structure for the material in that layer: 

<F (r) = J2 K (rRfco, . j = 1,2 . . . M (5) 

71 = 1 

where Nj is the number of bands that describe the band-edge bulk properties of the material in the j-th layer, and 
M is the number of layers. The symmetry of the material and the number of bands can be different in each region. 
The Nj component envelope function ^(r) = {^(r)}^! is slowly varying in the bulk-like region r € Rj, where it 
satisfies the k • p Schrodinger equation 

H j (r) = e* J '(r) , r e R j . (6) 

The Nj x Nj Hamiltonians £P are obtained after averaging of the perturbed microscopic Hamiltonian -ff m i cr (p + Tik) 
over a unit cell in each bulk-like regiono and includes terms up to the second order in the wave vector operator 
k = — iV. The boundary regions, however, are excluded from the averaging procedure, and the parameters of the 
resulting heterostructure MEMA Hamiltonians, H J , may have abrupt jumps from one layer to another. 

In what follows we study in details the boundary conditions imposed on the wave functions ^(rl = ^ (r) for r G Rj 
at the interfaces. The microscopic Hamiltonian H m i CT acting on <I>(r) is a self-adjoint operatoiO. This is equivalent 
to the conservation of charge in the full heterostructure region Q: 
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d 3 r|$| 2 =0. (7) 



If a/L is small, the contribution from the boundary regions can be neglected. Therefore, Eq. (Q), within the accuracy 
of the MEMA, is replaced by d/dt ("J, ^) = 0, where "3/ is the corresponding envelope function satisfying the BC and 
the inner product is now defined as 
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(8) 



This implies, in turn, that the heterostructure MEMA Hamiltonian, H = iP for r g Rj, acting on '5 has to be a 
self adjoint operator. Therefore for any two arbitrarily chosen functions ^ a and ^p, the condition 
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must be satisfied if the same BC is imposed on as well as on ^ p (see, for example, Ref. ^4|). The MEMA 
Hamiltonians H J , when defined for the ideally infinite homogenous bulk semiconductor, are self-adjoint, but when 
one adds the abrupt boundaries, the condition of Eq. (^|) may be violated by boundary terms. The requirement for 
H to be self-adjoint limits the choice of the boundary conditions that may occur, but there is an ample variety of 
BC that satisfy it. We shall see that any self-adjoint MEMA Hamiltonian preserves the normal component of the 
envelope flux density for an arbitrary state of the system. 

The MEMA Hamiltonians for the Nj component envelope wave functions have terms of first and second order in k: 
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where C 3 , each Bj^ (/i = x, y, z) and D 3 v (/i , v = x, y, z) are Hermitian Nj x Nj tensors of rank 0,1, and 2, respectively. 
They contain the energy band parameters for the material in Rj and are defined only in these bulk-like regions. Using 
Eq. ( |l0| ) we can write a velocity operator V = V J for r € Rj which acts on the heterostructure envelope functions: 



h dk 

A general envelope flux density matrix can be written: 



Ok 



J a/3 (r) 



(w Q ,V) 



r e R, , 



(11) 



(12) 



where "$> a and ^ p are two arbitrarily chosen solutions = ^ a (r 6 Rj), = &p(r 6 Rj) of the Schrodinger 
equation Eq. (||), defined in each bulk-like region, with energies s a and sp respectively. Noting that 



divJ Q/3 (r) 



(k* a ,V^)-(* Q ,kV^)- 



(V* a ,kf /3 ) + (kV* a ,f /3 ) 
and using Eqs. (fill nil) we arrive at the identity: 
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VcHVpj-lHVa,*/}) = -divJ^r), reR 



Thus, for the MEMA Hamiltonian H to be self-adjoint with BC imposed on \]/ a and \& p, Eq. (j|) requires 

Y^f d 3 r- divJ Q/3 = J2 I rfS J -J" /3 = 0, 



(13) 



(14) 



where dSj is an element of the surface bounding Rj. Noting that Eq. (|lj) holds for arbitrary surface shape, we find 
that local continuity of the normal envelope flux density matrix across each interface must hold: 



a/3 _ 



1 



= const . 



(15) 



Equation ( |15| ) is a generalization of tke-jaormal flux density conservation law, J T = J" a = const, that is often used 
as the starting point for deriving BCaO The condition J" 13 — const is a very general and strong requirement on 
the envelope functions at each point of the interface. For an unpcnetratable interface, or a semiconductor/vacuum 



4 



surface it reduces to — 0. We will consider this case only in Section IV, and shall assume, for now, that J" 13 7^ 
at the semicondutor heterointerface. 

The condition of Eq. (|l^) must be fulfilled independently of whether the number of the envelope function com- 
ponents on the two sides of interface is the same. For clarity in the following we will only consider heterostructures 
formed of semiconductors whose energy band structures are described by multiband Hamiltonians H J of same sym- 
metry and with TV,- = N in each region. For the more general case, a more complicated procedure similar to that 
suggested in Ref. QJJ is needed. The appropriate number of boundary conditions depends on the number of indepen- 
dent components of the envelope function ^> and normal envelope velocity V T — t ■ V^. If all the components of W 
and V T , are linearly independent on each side of the interface, the most general BC one can impose on the envelope 
functions are of the form: 

Xi ' ' T 2N ( * ) , (16) 



iV+ / ™ 1 iv 



where the indices "-" and "+" refer to neighboring points on neighboring bounding surfaces Sj and Sj+i, respectively, 
and T 2 n is the 2N x 2N transfer matrix. In the BC described by Eq. (|l6|) the material parameters of the MEMA 
Hamiltonian are included in the vectors of the normal envelope velocity. Equation (jig) is satisfied for all arbitrarily 
chosen a and (3 if and only if the BC of Eq. (Jig) are imposed on all \I/ Q as well as on all ^ p with the same T 2 n and 



TL V T 2N = " ^ , (17) 



->t / !n \ rp _ ( I N 
L ™ { -I N J l2N ~\-I N 

where In is the N x N unit matrix and the dagger denotes the Hermitian conjugate. Therefore, T 2 n must be state 
independent and characteristic only of the interface. For the conventional BC, ^ + = , = V~ , the transfer 
matrix reduces to the unit matrix, T 2 n = I 2 n- 

In the "pure" Kane models (D^ v = and H only contains terms linear in k) the N components of V T = f ■ B^P 
depend linearly on the N components of In this case, the general BC can be written 

= Tat*- , T N BtT N = B~ . (18) 

Here B T = f ■ B and Tjy is an N x N state independent transfer matrix. One can see that if ^ B~ , the transfer 
matrix Tjy can not be the unit matrix In- The number of independent equations in Eq. ( [l8|) is determined by the 
rank of the matrices B T and can be less then N. 

The elements of the state independent transfer matrix may depend on such properties of the heterojunction as band 
offsets and surface crystal symmetry due to, e.g., its reconstructions. Additional restrictions on the-,cpmponents of 
the transfer matrix can be obtained from symmetry consideration of the bulk and surface propertiesliilu in particular 
cases. For heterostructures which geometrical shape has some elements of symmetry one can expect the boundary 
parameters to be the same at the symmetrical points of the interface. For example, in a QW the parameters of the 
transfer matrix are the same along the plane, and in the spherically shaped heterostructure, the parameters of the 
transfer matrix are the same anywhere on the surface. 



III. ELIMINATION OF SPURIOUS SOLUTIONS FROM MULTIBAND K P MODEL: TRUNCATED 

BOUNDARY CONDITIONS 

Consider a heterostructure with a single heterointerface described by a piece- wise N x N Hamiltonian of Eq. ([h]), 
and with nonzero energy dispersion matrices D^. The envelope functions can be expanded in a complete set of N 
plane wave eigenf unctions, ^>i(ki), of the N x N MEMA Hamiltonian: 

JV 

* = ^CA(fc,), (19) 

i=l 

where the N values of kf are the roots of the bulk energy dispersion equation det[H(k) — e] = 0, and Ci are 
numerical coefficients. The set of 27V boundary condition equation in Eq. (|T7| ) can be solved only if the expansion 
in Eq. (|l9|), for both sides of the interface, includes all N of the plane wave eigenfunctions. If any, say m roots \ki\ 
(i = N — m + 1 , . . . N) lie outside the first Brillouin zone, the corresponding spurious solutions must be excluded 
from the wave function. In this case, 2m components of the envelope wave function and their velocities are linearly 
dependent on the other components, and only 2(N — m) independent BC can be imposed on the truncated function 
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$ = YliLi™ 1 Ci^i(ki). In this section we suggest an approximate procedure which allows us to derive a truncated set 
of general BC that must be imposed on the envelope function in order to satisfy J"^ = const. For the sake of clarity, 
we first consider a two band model (N = 2, m — 1) and a planar heterointerface. We then generalize it to arbitrary 
N and m and an arbitrarily shaped heterointerface. 

The Hamiltonian of the two band Kane model is written: 



H 



E c 
E„ 



H 



K 



H 



K 



2mo 



I 



\ 



a c kj 

~i~^k z 
n 



i~rk z 

n 

-a v kl 



(20) 



where k z = —idjdz is the wave vector along the normal to the interface, tuq is the free electron mass, P is the Kane 
interband matrix element, a c and a v describe the contribution of the remote bands to the electron, m c , and hole, ra„ 
effective masses respectively, and the distance between the bottom of conduction band, E c , and the top of the valence 
band, E v , is the energy gap E g = E c — E v . The band edge effective masses can be written as mo/m c ,„ = a c ,. u + E p / E g . 
Here E p = 2P 2 /tuq is the Kane energy which also characterizes the nonparabolicity of the electron and hole bulk 
energy spectra. Although this two band model completely neglects spin-orbit interactions, it describes the electron and 
light Tjinlfj band coupling in real semiconductors and is often used in one dimensional heterojunctionslill and spherical 
dotsEjOEj. For most semiconductors, \a c \ ~ |a„| ~ 1 (they can be either positive or negative) and E p w 20 ± 5eV. 
In the "pure" Kane model a c = a v — 0. 



For a planar heterointerface, the general normal flux density continuity condition, Eq. fll5|), takes the form: 



J. 



(z) = \[(*,V 9 ) 



(V*, <&)] = const . 



(21) 



where "J and $ are two arbitrary chosen eigenfunctions of the Hamiltonian in Eq. (p0|), = V z $> and V$ 
the respective normal velocities calculated with the two band Kane model velocity operator V z 



V z $ are 



m 



The gener al BC, Eqs. (H, 0), now reduce to 



.P 



-i— -a v k z 
n 



(22) 



iVJ, 



— T tr 



h 

-h o 



Ttr — 



h 

-h o 



(23) 



where T tr is an energy independent 4x4 transfer matrix, and r (I) refers to a point to the right (left) side of the planar 
interface. The conventional BC ( x I' r = ^ l and V£ — V$) are described by the unite transfer matrix T tr = I4. In the 
general case, a c a v 7^ 0, the BC of Eq. (^3|) represent 4 linear independent equations for two (conduction and valence 
band) wave function components ty v and two normal velocity components V* c , V$, v . In order to solve them, the 
envelope function of a state with energy e is written as a linear superposition of two plane wave eigenfunctions of the 



Hamiltonian, Eq.(EOJ): $ = Ci^'i(fci) + C^^^), where k\ 2 are the solutions of the bulk energy dispersion equation 
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Ei-e 
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2 J.2 



h 2 k 2 , 



(24) 



Generally we are only interested in states whose energy is either on the order of or smaller than the semiconductor 
energy gap. Using the condition E p /\a c a v \ ^> E g , \e—E c \, \e—E v \, we find for the two roots of the dispersion equation: 



h 2 k\ 



2{e-E v )(e-E c )m 
(e - E v ) a c + (E c - e) a v + E p 



n 2 k 2 



2E p mo 



(25) 



For typical values of the bulk energy band parameters, the second wave vector |fc 2 | ~ 1 A -1 lies beyond the first 
Brillouin zone in which the MEMA is valid. The corresponding eigenvector, ^2(^2) is a spurious solution (either 
evanescent or propagating, depending on the sign of the product {a c a v )EJi) and must be excluded from the expansion 
of the total wave function. The continuity of the normal flux density, Eq. (^l|), is now limited to slowly varying 
functions of ki, * = *i(fci) and $ = $i(fci), and leads to J??(z) = J* 1 $1 (z): 
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to 
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£$ — Ey 



E r 



Ey 



E, 



Ey — -Et, 
E„ 



P 
£<3> 



2rrioE p 



(fc?(e*) +fc?(e*) 



-Em 



a c a t 



2m E,, 



Const . 



(26) 



Here and e$ are the energies of the states described by the two component envelope functions ^(fci) = — 
^(fci), = ^(fci)) and $i(fci) = (<& c = 3>c(fci), $1; = ^(fci)) respectively. The two components of the normal 
velocities Vy c (V$ c ) and Vy v (V#„) are expressed through the wave function components ^ c ($ c ) and 'J,, (<&„). Since 
the envelope functions 'J and <£> are now truncated, the normal velocity and envelope function components are no 
longer linearly independent. Therefore, the Eq. (26) can not be fulfilled if a c a v ^ 0. (Fo r the case a c a v = the 
spurious solution does not arise, and is then the full solution of the Hamiltonian (20)). However we shall see 
below that we can impose an alternative set of the state independent BC that retain the continuity of the normal flux 
density J** for truncated wave functions within an accuracy consistent with the MEMA. We consider a J T ,i that 
differs from J Tj i of Eq. ( p6|) only by the terms 7 2 ~ a c a v which are small: j 2 = a c a v (£ — E c )(e — E v ) / E 2 w (fci/fe) 2 - 
A typical value of the momentum ki in a heterostructure with size L is fci ~ 1/L. Substituting the value of from 
Eq.(22) into 7 we find that j 2 ~ (ao/L) 2 <C 1. We can write then the approximate condition 
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E, 
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1 + a v 



E c — Ey 



= Const . 



which is fulfilled exactly if the envelope functions W and $ satisfy the truncated set of BC: 
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(27) 



(28) 



where the components of the 2x2 transfer matrix T tr are real and state independent and det[T tr ] = 1. One can 
see that the connection between the left and right hand side components of ^ c and ^ v has now energy dependent 
coefficients. This is a consequence of the linear dependence of the components of the truncated wave function and its 
normal velocity on the same side of the interface. However, the new energy independent transfer matrix T tr is still 
characteristic of the interface only. The truncated general BC of Eq. (|2|) no longer allow the simultaneous continuity 
of the both envelope function components ^(fci) and ^> v (ki) at the interface. For the case where the off diagonal 
matrix elements of T tr are negligibly small, ti2 ~ ?2i ~ 0, the general BC reduce to: 



P a 1 + a., 



E r 



E„ 



* c = Const . 
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ty v = Const , 



(29) 



(21 28) are satisfied exactly, 
with a c — a v — 0; also the 



which is similar to the BC of Eq. (Q) for heterostructures with parabolic bands. 

Equations (|2^, |2^) are general BC for the truncated smooth wave functions in semiconductor heterostructures with 
non-zero a c a v . These BC do not depend on the sign of a c a v (i. e., whether the spurious solution is evanescent or 
propagating). They are also valid for any Kane model with a c a v — 0, in which Eqs. 
The boundary condition Eq. ( p9| ) was first discussed in Ref. || for the " pure" Kane mode 
model with a c = 0, a v ^ was considered in Ref. 

To generalize this procedure and derive general BC for truncated slowly varying envelope functions for heterostruc- 
tures described by an N x N Hamiltonian having to spurious solutions, it is important to trace from where these 
solutions arise. Examining Eqs. (17,^5]) (and also see Ref. [II] which describes spurious solutions in the 8x8 model) 
we find that the spurious solutions \ki\ (i = N — rn + 1 , . . . N) of the dispersion equations for the Hamiltonian of Eq. 
( |l0| ) are proportional to the large components of the matrices B and inversely proportional to the square root of the 
product of the tensor D LIU conduction and valence band components. For the sake of definitenessEZI, let us assume that 
the matrices B contain the submatrices B = B 7 /7 of rank 2to, where 7 ~ <Xq/L is a dimensionless small parameter 
responsible for the spurious solutions. For example, in the two band model considered above the normal component 
BJ = f ■ B 7 is: 
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(30) 
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We can write the velocity operator of Eq. (11) as 



V = Ibt + V , (31) 

7 

and write the truncated slowly varying envelope function as the superposition of two linearly independent orthogonal 
wave functions \I/ = + ^ 2l such that the 2m component function ^ and n = N — 2m component function ^> 2 
satisfy the condition 

B^l^O, B^ 2 = 0. (32) 

Furthermore, we can write the function f ■ Vo\l/ = Vi + V% as superposition of a 2m component function V± and an 
n component function V% belonging to the same subspaces as the functions VPi and '$> 2 , respectively, i.e., (Wi, V 2 ) = 
( 1 t>2, Vi) = 0. Note that \& 2 is absent (n — 0) in the two band model; ^1 is the total slowly varying envelope function 
and V\ — Vo^i (Vo is the diagonal part of the velocity operator in Eq. ((H)). However, the functions V\ and V 2 do 
not coincide with f • V \l/i and r • V \p 2 , respectively, for the general case n^O. The flux density continuity of Eq. 



( |l5| ) can now be written as: 



J r " " = J* f + = Const , (33) 



T al3 1 ,~ , 1 



J^=-{^ la ,B^ ip ) + [(MJ la ,V^) + (V la ^ ip )} , (34) 
7 2 

Jr,2=\[^2 a ,V 2(3 ) + (V 2a ,^ 2 p)] . (35) 



2 

Neglecting terms of second order in 7, Eq. ( ^) can be expressed 



J^ = -[^ Xa ,B^ 10 ) , ^^^i + ^-Vi, (36) 



where V\ is, to this order, included in Equations ( |34j ) and Eq. (|36|) are the multiband extension of Eqs. ( |26| ) 
and ( ^7| ) obtained for the two band model. The vector {^1, ^ 2 , V 2 } has in total 2m + 2n = 2(N — m) independent 
components and the general boundary condition satisfying J"' 3 w J"' 3 = J"f + J" 2 = const can be written as: 



T 2m T12 



*2~ = ^2(JV-m) I *2 I ) T 2 (N-m) = ( ^ m ) • (37) 

zV+ / \ iV 2 ~ I 

We shall neglect the off diagonal matrices, i.e., set T\ 2 = T%\ = 0, (thus satisfying the stronger condition that both 
J"' 3 = const and = const) and obtain the following restrictions on the elements of the energy independent 
transfer matrices T 2m and T 2n : 



jit 1 D7 

2m 7 T 



T 2m = -B? 
s,+o 7 



Equations ( |37| ) and (|38|) represent 2{N — m) boundary condition equation for the truncated envelope wave function. 
These general BC do not allow the continuity of all envelope function components at the interface if BJ/j is not 
continues. While the elements of T 2m and T 2n are state independent, the function may be related to Vf^ by the 
energy dependent expression. ._. 

This procedure can be directly applied to heterostructures described by the 8x8 Pidgeon and Brown HamiltonianEa. 
The bulk energy dispersion equations for this Hamiltonian has 8 solutions, fc?, (taking spin degeneracy into account), 
for each energy e; two of these are spurious, and proportional to the large interband Kane matrix element P (see, for 



example, Rcf. 11). In the case of zero spin-orbit interaction, must include coupled electron and light hole band 
envelope functions, while ^> 2 includes the heavy hole band envelope functions only. For planar heterostructures Eqs. 
( p7| , |38| ) take into account the mixing of light and heavy holes for the states with finite in-plane momentum. However 
they neglect the effect of the " low interface symmetry"Ll that mixes light and heavy hole plane waves normally incident 
on the interface (because we assume that T 12 = T 21 = 0). The detailed application of our procedure to the 8x8 
models for planar and spherical heterostructures will be presented elsewhere. 
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IV. ONE DIMENSIONAL QUANTUM WELL WITH INFINITE POTENTIAL BARRIERS. 



Let us consider an effect of the general BC on the energy spectrum of a one dimensional quantum well having 
infinite potential barriers. The conventional BC for structures with an impenetrable barrier require vanishing of the 
envelope wave function at the barrier surface, but this requirement has never been justified (see for example Ref. 
|39| ). In general, the self-adjointness of the MEMA Hamiltonian requires vanishing of the normal flux density at the 
boundaries: 

J?P(z = ±LtO)=0, (39) 

where 2L is the thickness of the quantum well. If we assume that the two interfaces are completely identical, the 
quantum well possesses reflection symmetry about z = and all solutions of the Schrodinger equation are characterized 
by their parity according to whether their components are even or odd under reflection. However, the BC are local 



characteristics of each interface. Therefore, the Eq. (39) must be fulfilled at each interface independent of the 
symmetry of the functions ^ a and ^p. 

Now we consider the effect of the general BC for the case of the two band Kane model. In order to satisfy Eq. ( |39| ) 
in a quantum well with symmetric interfaces, one can write the general BC for the conduction ty c and valence band 
tyy components of the envelope wave function (which do not contain spurious solutions) as 

* C (±L T 0) (l - «« £ ^ £ ) = TO * V (±L T 0) (l + "c 1 ^) > (40) 

where is a real number and the sign difference is due to the opposite parity of the conduction and valence band 
components of the envelope functions. The surface parameter Q does not depend on the energy or symmetry of either 



state. Equation (|4Cj) can be derived directly from the general BC of Eq. (|28|) obtained in Sec. Ill for the case of 
a finite potential barrier. To do this one has to assume that ^ = outside the quantum well and use the state 
independent transfer matrices T t +, for z = L — 0, and T t ~, for z = —L + interfaces, which differ only by the sign 
(opposite) of the offdiagonal matrix elements and satisfy the condition det[T t +] = det[T t ~] = 0. 

Neglecting j 2 ~ (clq/L) 2 terms, Eq. ( |Io| ) can be written as two separate equations for the conduction and valence 
band components, respectively 

$ c (±LT0)=±T c a ^-<(±LT0), T c a = 6 A (41) 
m c (E e ) 2P 

^ v (±L T 0)=±T v a -^-%(±L T 0) ) 7> = i^. (42) 
m v (h h ) 6 2P 

Here the energy-dependent effective masses in the conduction and valence bands are given by 

mo ---l. m " -a v + , (43) 



m c {E e ) c E e + E g ' m v (E h ) v E h + E g ' 

where the electron, E e , and hole Eh energies are measured from the bottom of conduction and the top of the valence 
bands, respectively: E e = e — E c and Eh — E v — e. The even, (+), and odd, (-), solutions to Eqs. (38,39) can be 
written 

*c(«) (*) = Mv) cos(# (t)) Z), K(v)( Z ) = A c(v) sin (C(„) z ) . ( 44 ) 

where A CiV are the normalization constants. We can derive equations for the energies of the even and odd electron 
and hole quantum size levels: 

2m L 2 m c (Et) \ E p ) h 2m L 2 m v (E^) \ ° E p J V ' 

where </>^ and are the solutions of the equations: 

^[tanfefrF 1 = T—^—^tJ— > ^[tM^tr 1 = T !^l-i- . (46) 
mo J c ao mo Ivao 

Figure |l| shows the dependence of the two lowest quantum size electron levels, E£ and E~ , on the well width L 
as a function of of the surface parameter T c a$. In these calculations we use for m c the electron effective mass at the 







bottom of the conduction band m c — O.lmo, and a band gap energy E g — 1.7eV which are close to the parameters of 
CdSe and CdTe. We compare the effect of the general BC on the quantum size levels in the parabolic EMA (Fig. la), 
in the "pure" Kane model with a c = a v = (Fig. lb), and in the "full" two band model with a c a v ^ (Fig.lc). 

Equations ( fl5| , [l6|) describe the energy of the quantum size levels for coupled conduction and valence bands. However, 
the conduction and valence band energy spectra can be considered as being independent when the energies of the 
electron or hole levels are much less than the energy gap, |-E e |,|-E/j| <C E g ,E p . This limit case is realized in thick 
quantum wells and is described by the equations for the simple parabolic bands that one obtains by neglecting the 
energy dependence of the effective masses and the E ey h/E p terms in Eqs. ( |45|j46| ). For this case the surface parameters 
T c ao and T v ao can be chosen independently, and the conventional BC ^f c (±L) = are realized for T c ao = 0. Fig. 
|l|(a) shows the effect of T c clq ^ on the electron quantum size levels in the parabolic EMA. One sees that positive 
and negative values of T c ao shift the quantum size levels energy up and down, respectively, from those obtained using 
conventional BC. The effect is negligible when \T c ao\/L <C <j>^r mo/m c but become noticeable in narrow wells and is 
greater for higher energy levels. 

Fig. 0(b) shows the effect of the general BC on the electron quantum size levels in the pure Kane model. One sees 
that the effect of T c ao ^ is very similar to that in the parabolic EMA. The size dependence differs only in the nature 
of the nonlinear dependence on 1/L 2 . This is because the Kane model takes the nonparabolicity of the conduction 
band into account. 

The size dependence of the electron levels calculated with T c ao = is shown in Fig. |](b) only for comparison. 
The conventional BC do not hold in the Kane model, because the surface parameters for coupled conduction and 
valence bands are related by (T c ao)(T v ao) = h 2 /2E p m a . Choosing T c a = for determining the electron energy levels 
corresponds to choosing T v a — > oo, which does not describe the hole energy levels. The condition |T c a | = |T u a | = a* 
realized for |0| = 1, describes a symmetric (relative to the to the center of the band gap) structure of the electron 
and hole quantum size levels in semiconductors with a c = a v . Although these " symmetric semiconductor structures" 
do not exist in nature, the parameter a* = h/2P w 0.45 ± 0.06 A gives a reasonable value of |T c ao| and |T„ao| in real 
semiconductor structures. If symmetric BC hold 

i too a* i TOn a* 

«1, 4>t ° — <1, (47) 



m c (E e ) L v m v {Eh) L 

and the solutions of Eq. ^6| for the lowest electron and hole quantum size levels, 4>^r v , are close to those given by 
T CjV ao — 0. 

Figure (c) shows the size dependence of the electron levels calculated in the full two band model with the symmetric 
surface parameter \Tao\ = a* = 0.43 A. One can see that the term linear in a v in Eq. pq ) for the electron energy 
levels becomes important when they comparable with the band gap energy. On the other hand changing the sign of 
a v leads only to small changes of the level energy. 

The dependence of the lowest electron quantum size levels on the surface parameter T c ag for the Kane and full two 
band models for wells with L = 30 A and L = 12 A, respectively, is shown in Fig. ^. The surface parameter is varied 
from — 3.0A to 3.0A for L = 30 A, and from — 1.2 A to 1.2A for L = 12A, respectively, so that |T c ao| < Lm c /mo 
is fulfilled. One can see that varying the surface parameter in this range produces a monotonic change of the first 
even and odd electron levels. The difference between models with different a c and a v (for the same m c ) is small for 
L = 30 A, but can be important for L = 12 A. 

It is interesting to note that for positive values of the surface parameter O (T c ao > 0, T v ciq > 0), equations (45 4^ 



may have even and odd solutions with an energy in the forbidden gap: E e < 0, Eh < 0. In wide wells the energy of 

'+ — F~ 

(e,v),S ~ (e,v),S 



these gap levels do not depend on their symmetry, Ef v \ s = E, v \ s — Er etV \g, and can be found from 



h 2 m c (E e:S ) 1 / E St s 



2to to (T c a y \ E p 



h m v {E h . s ) 1 / E h .s \ 

E h ,s = -7, 7^ — >2 • U-Q^-^r- • (49) 

2to to (T v a y \ E p J 

In wide "symmetric wells" these gap states have the same energy E e ^$ = Eh.s = ~E g /2, and are localized within 
a layer of thickness a s — 2a* E p / E g near the surface. These solutions contain no contributions from the unphysical 
spurious solutions, and thus are not artifacts of the k • p model. 

V. DISCUSSION AND COMPARISON WITH OTHER APPROACHES 

The occurrence of discontinuities in the envelope wave functions at the heterointerface is one of the most important 
consequence of general BC. These discontinuities have a strong effect on the mathematical procedures often used for 
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the calculation of various physical properties of heterostructures having finite potential barriers. In these procedures, 
the Schodinger equation is integrated across the interface or a Fourier transformation is performed using the piece- 
wise spatially determined material parameters with the help of generalized step functions. However, the integration 
of terms such as products of ^k z or Pk z with envelope wave functions, 'J, which are discontinuous across an interface 
may lead to mathematical uncertainties: integration, .of the product of a step function and a (^-function is not well 
defined. To resolve this problem, a non-uniqueaE9li§0 symmetrized form of the kinetic energy operator is-used, 
and the BC for the envelope function 4" are obtained by requiring that Hk^> be integrable across the interfacdfErEil. 

For exampje, in order to obtain the general BC of Eq.(4) for parabolic bands, symmetrized kinetic energy operator 
of the formEl 

= < 50 > 



is used. The case a = then leads to the conventional BC of Eq. (j^l^Ihe same symmetrized form, k z — k z , is usually 
used for the diagonal terms of the multiband kinetic energy operatoroEj. To derive the BC of Eq. ( p9| ) for the "pure" 
two band Kane model (a c — a v — 0) that allow to integrate Schrodinger equation across the interface, one writes the 
kinetic energy operator as: 

~ mo ^ -P l -*k z P a ) ■ (51) 

Here, now, there is no value of a that gives the symmetriKed-form l/2(fc z P + Pk z ), that is usually used for the off 
diagonal terms (linear in k) of the multiband Hamiltonianfi-3'EJ An asymmetric ordering, corresponding to a = 1 was 
suggested in Ref. |33| for the model with a c — 0. The a dependence of Eqs. (|5^) and (pi"| ) clearly demonstrates that 
it is the BC that determine the integration across the abrupt heterointerface and not the other way round. 

Fourier transforming the Schrodinger equation has been suggested by Winkler and RosslerEJ as an alternative ap- 
proach to the MEMA problem for heterostructures with finite potential barriers. Then one does not deal explicitly 
with BC when solving the resulting integral equations for the momentum space envelope function, and avoids un- 
physical spurious solutions by restricting the range of integration to |fe| <C 2ir/ao. However, the resulting form of the 
momentum space MEMA Hamiltonianj-depends on the particular symmetrization procedure chosen for the kinetic 
energy operator in configuration spaceEl, and, therefore, is again determined by the BC imposed on the envelope 
function. To illustrate this, we derive the explicit form of the momentum space Schrodinger equation for the one band 
EMA and the two band "pure" Kane model, using general BC. Fourier integration of the EMA Schrodinger equation 
with the kinetic energy operator of Eq. ( |50| ) leads to the following integral equation for the momentum space envelope 
function F(k): 

h 2 f°° ii i ii i ii 

— \ qk ' M- a {k - q)M 1+2a (q - k')M- a (k' - q)F(q)dq'dk dq + 

* J -oo 

/oo 
F(q)V c (k - q)dq = eF{k) , (52) 
-oo 

where F(k) satisfies the normalization condition \F(k)\ 2 dk — 1, and 

1 r°° i r 00 i 

V c (k) = — E c {z)exp{-ikz)dz, M u {k) = — / — — exp(-ikz)dz . (53) 
2tt J_ oa 2tt J_ 00 m v (z) 

Here E c (z) and m(z) are the energy of the bottom of the conduction band and the electron effective mass in each 
region of the one-dimensional heterostructure considered. The configuration space envelope function, in the j— th 
region of the heterostructure, is now given by 

1 r°° 

f 3 {z) = m- a -—= F(q)M- a (k-q)exp(ikz)dkdq, (54) 



'2tt J- 



and satisfies the general BC of Eq. (0). For the "pure" Kane model, Fourier transforming the Schrodinger equation, 
using the kinetic energy operator of Eq. (pl|), leads to the following coupled integral equations for the two component 
momentum space function {F c (k), F v (k)} : 
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/ F c (q)V c {k - q)dq + i— q'G a (k - q')G x ^ a {q' - q)F v (q)dqdq' = eF c {k) , (55) 

J-oo m J-oo 

-l — q'G 1 - a (k-q')G a (q'-q)F c (q)dqdq'+ F v {q)V v {k ~ q)dq = eF v (k) , (56) 



^0 J — og ^ — oo 

where (|F c (fc)| 2 + \F v (k)\ 2 ) dk = 1, and 

i r 00 i r 00 

V cv (k) = — E cv (z)exp(-ikz)dz, G u (k) = — P v (z) exp(-ikz)dz . (57) 

2?r J-oo ' 27r J-oo 

Here E v (z) and P(z) are the energy of the top of the valence band and the Kane matrix element in each region, 
respectively. The configuration space envelope function $ = {^ c ,^ v } in the j—th region of the one-dimensional 
heterostructure is now given by 

1 [°° 

W c (z) = Pf a ■ -= / F c (q)G a {k-q)exp(ikz)dkdq, (58) 



/2tt 

1 

V2tt J- 



*l{z) = P"' 1 ■ I F v {q)G!- a {k - q) exp{ikz)dkdq , (59) 



and will satisfy the general BC of Eq. ( p9| ) with a c = a v = 0. 

General BC, for heterostructures described by the MEMAJJamiltonian, that conserve the normal component of the 
envelope flux density J T have been suggested by Kisin et alER We have shown that the general requirement that the 
heterostructure MEMA Hamiltonian be self-adjoint leads to the more general condition J"" = const. That is, that 
J T must be continuous for arbitrary chosen envelope functions ty a , *S>p. This generalization to a ^ (3 is important, 
because it shows that the components of the transfer matrix connecting the components of the envelope function 
^ and envelope normal velocity V T on the two sides of the interface, are state independent and are characteristic 
properties of the interface, notj the states. This is not true for a transfer matrix which connects the components of 

and its derivatives $ (see Ref. fLl|), because the relationship between the components of V T and '5, \& , on the 
same side of the interface, in the multiJband model may depend on the symmetry (the value of the interface parallel 
momentum in planar heterostructuresEiJ and the total orbital angular momentum and parity in spherical dotsLj) or 
on the energy of the state. 

An explicit treatment of "interface" effects can ba. seen an several k ■ p models of heterostructures that consider 
exact electron wave functions in the interface region&r H a. Some of these models result in equations in addition 
to the bulk k • p equations that describe the interfaceaa. The parameters of , these " interface" Hamiltonians must 
be determined from microscopic wave functions or obtained from the comparison with experiment. However, these 
advanced models are generally much more complicated than the MEMA, completed by an appropriate choice of 
boundary conditions, which has,-already successfully described such interface effects as the T — X intervalley mixing 
in GaAsLAlAs heterostructuresQ'El, the heavy-light hole plane wave mixing at zinc blende interfaces under normal 
incidenceQ and the effect on the interband light absorption of the intervalley conversion of the electron at the surface 
of indirect-band-gap semiconductors^. 

Heterostructures with abrupt heterointcfaces, (L ^> a), are described by the MEMA models using energy inde- 
pendent transfer matrices which characterize the effects of the interfaces on the carrier wave functions within some 
energy interval. The energy of the state considered should be within the range of validity of the chosen bulk MEMA 
Hamiltonian in each layer. This means, that the energy should be smaller than the distance to the energy band ex- 
trema not explicitly included in the multiband model. In the appropriate energy interval, the elements of the transfer 
matrices can be treated as energy independent trial parameters in fitting the experimental data. 

The number of independent transfer matrices required to ensure the flux continuity of Eq. ( |i"5| ) , is determined by 
the number of physically nonequivalent interfaces of the heterostructure. Therefore, to write general BC, one must 
model the symmetry of the interface, i.e. the interface geometry and the symmetry of the material properties on the 
interface. On the other hand, this is true for any explicit determination of the full microscopic wave function near 
the interface. 

In conclusion, we have derived a general form of state-independent BC for multicomponent envelope functions that 
are valid under the same conditions as the MEMA itself. Spurious components of the wave functions are eliminated 
by requiring that the envelope flux density be determined to the same order of approximation as the MEMA itself. 
The effect of the general BC on the electron and hole energy spectra has been demonstrated in the two band model 
for a potential well with infinite potential barriers. This procedure, using the general BC for truncated envelope wave 
functions, can be applied to any MEMA model containing any type of spurious solution and to heterostructures with 
finite or infinite potential barriers. 
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FIG. 1. The size dependence of the two lowest, even (E^) and odd (E~), electron quantum size levels in a well with infinite 
potential barriers calculated as a function of the surface parameter Tao: (a) in the parabolic EM A model; (b) in the pure 
Kane model (a c — a v — 0); (c) in the full two band model (a c a v 7^ 0). In (a) and (b), the solid, dashed and dotted curves 
correspond to Tao = 0, Tao = 0.5 and Tao = —0.5 A, respectively. In (c) the solid and dotted curves are calculated with 
a v — —2 and a v = 2, respectively, upper curves for each level correspond to Tao = 0.43 A and lower curves to Tao = —0.43 A. 
Other parameters used in calculations are shown in the Figure. 
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FIG. 2. The dependence of the two lowest , even (E^) and odd (E~), electron quantum size levels on the surface parameter, 
T c ao, calculated in a well of width 2L surrounded by infinite potential barriers: (a) L = 30 A and (b) L = 12 A. The solid 
line is for the pure Kane model with a c = a v = and E p — 17.0 eV, the dashed line is for the symmetric two band model 
with a c — a v = —2 and E p — 20.4 eV, and the dotted line for the asymmetric two band model with a c = —2, a v = 2 and 
E p — 20.4 eV. Other parameters used in the calculations are shown in the Figure. 
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